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Abstract. In the absence of inelastic scattering, Feshbach resonances produce poles 
in scattering lengths and very large peaks in elastic cross sections. However, inelastic 
scattering removes the poles. Whenever the resonant state is coupled comparably to 
the elastic and inelastic channels, the scattering length exhibits only a small oscillation 
and peaks in cross sections are significantly suppressed. A resonant scattering length 
is defined to characterize the amplitude of the oscillation, and is shown to be small 
for many collisions of ultracold molecules. The results suggest that cross sections for 
some ultracold collision processes will be much less sensitive to details of the potential 
than has been expected. 
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1. Introduction 

A Feshbach resonance pQ occurs when a bound state of a 2-particle system lies above 
a dissociation threshold and is coupled to the continuum. Collision properties show 
sharp features (peaks and troughs) near the energy of the resonance. In recent years, 
Feshbach resonances have come into prominence in the study of ultracold atomic gases. 
In these systems the positions of resonances can often be adjusted using applied magnetic 
fields, and it is possible to control the interactions between atoms by tuning resonances 
to near-zero collision energy [2j [31 0]. Magnetic tuning through Feshbach resonances 
has been used to produce molecules in both bosonic and fermionic quantum gases. 
Long-lived molecular Bose-Einstein condensates of fermion dimers have been produced, 
and the first signatures of ultracold triatomic and tetraatomic molecules have been 
observed. The new capabilities in atomic physics have had important applications 
in other areas: for example, the tunability of atomic interactions has been used to 
explore the crossover between Bose-Einstein condensation (BEC) and Bardeen-Cooper- 
Schrieffer (BCS) behaviour in dilute gases. There is now great interest in extending 
the capabilities from ultracold atomic to molecular systems, to explore the properties 
of dipolar quantum gases and develop new forms of quantum control. 

Most interpretations of Feshbach resonances have used concepts from the two- 
channel model [I], in which the bound state and the continuum are each represented by 
one scattering channel. This captures much of the crucial resonant behaviour observed 
in ultracold atom-atom scattering. In particular, it predicts that the scattering length 
passes through a pole and the elastic scattering cross section exhibits a very large peak 
at a zero-energy resonance. However, it is known from early work on nuclear reactions 
[5] that inelastic processes suppress resonant peaks in cross sections. The purpose of this 
paper is to explore the consequences of such effects for ultracold atomic and molecular 
collisions. Whenever the resonant state is coupled comparably to the incoming and 
inelastic channels, the scattering length exhibits only a small oscillation and the peaks 
in cross sections are dramatically suppressed. This is particularly important for the 
prospect of controlling molecular collisions. 

This paper will first summarize the results of 2-channel resonance theory, to define 
notation and establish a basis for comparison. The major differences introduced by 
inelastic scattering will then be considered. The results are general, but to assist 
visualisation the equations will be illustrated with examples taken from the elastic and 
inelastic scattering of NH molecules with He [6] . 

2. Resonances in the absence of inelastic scattering 

When there is only a single open channel with orbital angular momentum I, the long- 
range wavefunction may be written 




(1) 
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where 5(k) is the phase shift and the wave vector k is defined in terms of the kinetic 
energy E kin and reduced mass [i by E^m = h 2 k 2 /2fi. In the ultracold regime, cross 
sections are dominated by s-wave scattering, with I = 0. The most important parameter 
is the energy-dependent s-wave scattering length a(k), defined by 

a(k) = " tan / (fc) . (2) 

This becomes constant at limitingly low energy, with corrections given by effective range 
theory [TJ, 

o(jfc) = a(0) + h 2 r a(0) 2 + 0(k 4 ), (3) 

where r is the effective range. The elastic cross section is given exactly in terms of a(k) 
by 

For collisions of identical bosons, the factor of 4 is replaced by 8. However, the present 
work will omit such extra factors of 2. 

If there is only one open channel, the behaviour of the phase shift 5 is sufficient to 
characterize a resonance. It follows a Breit-Wigner form as a function of energy, 



5(E) = dbg + tan 



(5) 



_2(E rcs — E) 

where 5^ g is a slowly varying background term, E res is the resonance position and T B 
is its width (in energy space). The phase shift thus increases sharply by n across the 
width of the resonance. In general the parameters 5b g , E Tes and Te are weak functions 
of energy, but this is neglected in the present work apart from threshold behaviour. 

As a function of magnetic field at constant E k i n , the phase shift follows a form 
similar to Eq. [5], 



5(B) = 5b g + tan 



-i 



(6) 



_2(B TCS — B) 

where B res is the field at which E Tes = E = -E t hresh + -^kin- The width T B is a signed 
quantity given by Tb = T^/A/i, where the magnetic moment difference A/x is the rate 
at which the energy -Ethresh of the open-channel threshold tunes with respect to the 
resonance energy, 

A ^-^thresh G^res , > 

A »=^IB -IB' (7) 

T B is thus negative if the bound state tunes upwards through the energy of interest. 

Across an elastic scattering resonance, the S matrix element S = e 2lS describes 
a circle of radius 1 in the complex plane as a function of either energy or magnetic 
field, as shown in the left panel of Figure HJ In the ultracold regime, the background 
phase shift <5b g goes to zero as k — >• according to Eq. [2] (with ab g constant and finite), 
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Figure 1. The resonant circles described by S matrix elements for low-energy elastic 
scattering for two different resonances in He + NH ( 3 X~). Left panel: the circle of 
radius i when only elastic scattering is allowed: incoming channel n — 0, m s = — 1 at 
Eu-a = 10~ 6 K. Right panel: the much smaller circles (note the different scale) when 
both clastic and inelastic scattering are allowed: incoming channel n — 0, m s = at 
Eu-a — 10 -6 K (green, smaller circle) and 4 x 10~ 6 K (red, larger circle). The crosses 
show values far from resonance. In both cases the resonant state has n = 0, m s = +1. 



but the resonant term still exists. The scattering length passes through a pole when 
5 = (n + |) 7T, corresponding to S = — 1. The scattering length follows the formula [8], 

a(B) = a hg [l - Ab 1 . (8) 

The elastic cross section given by Eq. H] thus shows a sharp peak of height An/k 2 at 
resonance. The two widths r# and A# are related by 

T B = -2a hg kA B . (9) 

At limitingly low energy, T B is proportional to k [2] while A B is constant. 

3. Resonances in the presence of inelastic scattering 

In the presence of inelastic collisions, the scattering matrix has elements <S#/. The 
diagonal S-matrix element in the incoming channel has magnitude Soo < 1 and may 
be written in terms of a complex phase shift 5q with a positive imaginary part [9J, 

Soo(ko) = e 2i5o(fco) , (10) 

where k is the wave vector in the incoming channel. This can be expressed in terms of 
a complex energy- dependent scattering length, a(ko) = a{ko) — if3{ko) [HI [11], defined 
by analogy with Eq. [2] as 

-tan5 (fc ) 1 fl-S 00 (k )\ 
k ik \1 + S 00 {k ) J 
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a(k ) again becomes constant at limitingly low energy. The elastic and total inelastic 
cross sections are exactly [12] 



47T 


a 


2 


1 + *3I 


a 


2 





and 



° fco) " ' (13) 

When there are several open channels, the quantity that follows the Breit-Wigner 
form ([5]) or ([B]) is the S-matrix eigenphase sum [T3J [Hj, which is the sum of phases 
of the eigenvalues of the S matrix. The eigenphases and the eigenphase sum are real, 
unlike the phases <5j obtained from individual diagonal elements, because the S matrix 
is unitary, so that all its eigenvalues have modulus 1. 

Across a resonance, the individual S matrix elements describe circles in the complex 
plane [13 [16], 

Sw(E) = S hg , u > - £_|^ilV2' (14) 

where g B % is complex. The radius of the circle in Sa' is \gEigEi>\/^E- The partial width 
for channel i is usually defined as a real quantity, Y Ei = \gEi\ 2 , but here we also need 
a corresponding phase fa to describe the direction of the circle in the complex plane, 
g Ei = TEie 21 ^. For a narrow resonance, the total width is just the sum of the partial 
widths, 

T E = Y, T Ei- (15) 
As a function of magnetic field at constant -E^m, 

SAB) = S hs , u > - b _b2 9 +It b /2> (16) 

where gs% = gEi/^l^ 1 ^ 2 and the width T B and partial widths Y Bi are signed quantities, 
T B = r B /A/i and T Bi = r K /A/i. 

The partial widths for elastic channels (degenerate with the incoming channel) are 
proportional to k at low energy. We may define a reduced partial width 7^0 or 7^0 for 
the incoming channel by 

^Eo(ko) = 2fc 7i?o or T B0 (k ) = 2/c 7bo, (17) 

and the reduced widths are independent of ko at low energy. By contrast, the partial 
widths for inelastic channels depend on open-channel wavefunctions with large wave 
vectors ki and are effectively independent of k$ in the ultracold regime. If the inelastic 
partial widths Y B i (or T Bi ) are non-zero, they eventually dominate T B o (or T^o) as ko 
decreases. The radius of the circle (TT6|) described by Soo thus drops linearly to zero as 
ko decreases, as shown in the right panel of Figure [TJ This is qualitatively different from 
the behaviour in the absence of inelastic channels. 
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As a function of magnetic field, the scattering length passes through a pole only 
if 5 passes through (n + |) 7T, corresponding to Soo = —1- If there is any inelastic 
scattering, (T^ol < |Tb| and this does not occur. When the circle in £00 is small, the 
phase shift £0 and the scattering length a show only small peaks or oscillations across a 
resonance. 

The expression (fT2l) for the elastic scattering cross section saturates at a value 



(id « 4ir/k 2 when \a\ 3> k . Such values of \a\ occur only when 



8, 



-(n+l 



7T 



< 1 



and thus when T B is strongly dominated by T B0 . Since T B0 is proportional to k and the 
inelastic contributions T B i are independent of ko, there is a lower bound on the value of 
krj at which this occurs. Denoting the sum of inelastic contributions to T B as r£ el , this 
is given by 



-unel I 



< |r B0 | = 2/c |7bo| 



k 



> 



"pinel 
1 B 



27. 



BO 



(19) 



The radius of the circle in S o is T B0 /T B . For small k , where Eq. [IT] applies, this 
is approximately 2ko , y B o/r i B el . The formula followed by the complex scattering length 



is 



a(B) = a hg + 



2(B-B res )/Y^ + i 



(20) 



where resonant scattering length that characterises the strength of the resonance, 

Ao+fcoQbg) ^21) 



2 7B0 g 2i( 



res pinel 
1 B 

Both a res and the background term db g can in general be complex and are independent 
of k at low energy. The phase correction +2/c ab g in Eq- ED is needed to keep the phase 
of a res independent of ko. The explicit expressions for the real and imaginary parts of 
a(B) are 



a(B) = a hg + 



Or 



2(£? - A-esViT 1 +A 



(3{B) = p hg + 



2(B 

+ /9r, 



5 res )/r 



inel 
B 



2{B - B v 



+ 1 

0/r 



inel 
B 



"'inel 



(22) 



(23) 



2(b - s res )/n 

where a(-B) = a(-B) — i/9(-B) and similarly for a res and a\> g . The peak profiles for the 
elastic and total inelastic cross sections are given by Eqs. [12] and [T3j 

In the special case where the background scattering is elastic (ab g is real), unitarity 
requires that the circle in Soo must loop towards the origin. This requires that a res is also 
real. Across the width of the resonance, the real part a(B) of the scattering length a(B) 



then oscillates about db g by ±a res /2 and the imaginary part peaks at (3(B) 



When 



the background scattering is inelastic, however, be complex and the circle in Soo 

does not point directly towards the origin. The lineshapes are then unsymmetrical, and 
(3(B) (and hence the inelastic rate) can show a trough as well as a peak. Nevertheless, 
the overall magnitude of the oscillations in the scattering length is still governed by a res . 
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The behaviour derived here is analogous to that observed when laser light is used to 
tune scattering lengths [T~7] [TO] . However, in that case the amplitude of the oscillation 
depends on the ratio of excitation and spontaneous emission rates, which both depend 
on the same dipole strength (though the ratio of rates can be tuned with laser intensity). 
In the present case a ves depends on independent elastic and inelastic couplings. If a res 
is small, the resonant oscillations in cross sections and the scattering length are small. 

The results (1201) to (1231) are valid when fc 1 a res | <C 1. Whenever k a Tes j£> 1, Eq. [8] 
fails at values of \a\ small enough to reduce the height of the peak in the elastic cross 
section given by Eq. [T2j Conversely, when fc |a r es| 3> 1, «Soo describes a circle of radius 
close to 1 in the complex plane; the behaviour of the scattering length is then well 
described by a 2-channel model and the peak in the elastic cross section is of height 
~ Air/k^. 

The elastic partial width T^o is proportional to ko at low energy but becomes 
constant at high energy. It may be written [18] 

T B0 (k )=T B0 C (k )- 2 1 (24) 

where T B0 is independent of k and depends on the short-range coupling between the 
bound state and the incoming channel. The factor Co(fco) 2 is the amplitude matching 
function of multichannel quantum defect theory, which is 1 at high energy but near 
threshold is [18] 

C (k )- 2 = ha [l + (1 - a bg /a) 2 ] , (25) 

where a is the mean scattering length [TO], a = 0.478(2/iC 6 //i 2 ) 1/4 for a Van der Waals 
potential —Cq/t 6 . The transition between the linear and constant regimes depends on 
Cq and the reduced mass [20J, but typically occurs around E kin /k B = 1 mK. 

The height of the peak (or size of the oscillation) in the total inelastic cross section 
is proportional to |a res |. This in turn depends principally on the ratio of T B0 and Tg 61 . 
Two very different cases may be distinguished. If the same coupling term connects the 
bound state to the incoming and inelastic channels, it is likely that T BQ and Tg el will 
be comparable. Under these circumstances a rcs will be of the order of a and there will 
be relatively small oscillations in the scattering length. Conversely, if coupling to the 
inelastic (exoergic) channels is much weaker than coupling to the elastic channel, a res 
will be large and the scattering length will exhibit a large oscillation resembling a pole. 

It is important to realize that a res (and thus the strength of the resonance) depends 
on the relative magnitudes of the couplings from the resonant state to the elastic 
and inelastic channels. This is not necessarily the same as saying that the degree of 
suppression depends on the strength of inelastic scattering. 

The peaks in individual inelastic cross sections can be rather larger than those in 
of°ei, because the radius of the circle in S 0i is (2k \a ms \T Ei/T™ 1 ) 1 / 2 , which is considerably 
larger than 2/c |a res | for small ko- 
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4. Examples from low-energy atomic and molecular scattering 

For atomic collisions, the couplings to inelastic channels are sometimes weak enough 
that a 2-channel model remains accurate even when inelastic scattering is energetically 
allowed. For example, Donley et al. [21] and Thompson et al. [22] have produced 
85 Rb2 molecules by magnetic tuning in the vicinity of a Feshbach resonance between 
(f,rrif) = (2, —2) states of 85 Rb near 155 G. The (2, —2) state is not the lowest in a 
magnetic field, and the molecules can decay by spontaneous spin relaxation to atomic 
levels with f — 2 and m./ > —2. The resonant state has Mp = rrifi + m/2 — ~ 4, so 
this decay requires a change in Mp and involves very weak magnetic dipole coupling. 
However, the coupling between the resonant state and the incoming channel (also 
Mp = —4) is through much stronger central terms in the potential. Kohler et al. 
[23] have used coupled channel calculations including spin relaxation to characterize the 
resonance and obtained a^ g = —484.1 a and A# = 10.65 G. Their lifetime r = 32 /j,s 
for the bare resonance state corresponds to Tg el = h/rAp = 0.090 G. With these 
parameters, a res = 1.14 x 10 5 a . The temperature in the experiments of Thompson 
et al. [22] is 30 nK, corresponding to fco = 4.3 x 10 -4 <2g . In this system, therefore, 
&o a res ~ 50 and the resonant behaviour of the scattering length and the elastic cross 
section is well approximated by a 2-channel model. 

The situation is very different for rotationally inelastic molecular scattering, where 
the potential anisotropy couples the resonant bound state to both the incoming and 
inelastic channels. Under these circumstances a res will generally be small. In separate 
work, we have described numerical tests of the equations derived here for He + NH( 3 £~) 
scattering in a magnetic field [6j. This is a very weakly coupled system, and for the 
rotational ground state (n = 0) of NH the channels with different spin projections m s 
are coupled only indirectly via excited rotational levels. The background scattering is 
essentially elastic, so a^ g and a res are real. Fig. [2] shows the real and imaginary parts 
of the scattering length for magnetic tuning across an inelastic scattering Feshbach 
resonance in this system. Even He + NH, where the inelastic couplings are much weaker 
than in most other molecular systems, a res ~ 9 A and koa ves <C 1. The oscillations in 
scattering lengths and elastic cross sections are strongly suppressed at low energies. 

There are also atomic systems where the coupling to inelastic channels is strong 
enough to suppress the oscillations in scattering lengths. Such effects have been 
observed, for example, in calculations on collisions of Sr ( 2 P2) atoms [24j, where the 
bound state is coupled to both the incoming and inelastic channels by anisotropic 
potential terms. 

Eqs. ( 1201) to (1231 can be adapted to apply to any parameter A that tunes scattering 
resonances across a threshold. The ratio Yxq/Y™ 1 is the same for any such parameter 
(and is equal to Tpo/V™ 1 ). The resonant scattering length therefore has the same 
value for any parameter A. universal measure of the strength of a low-energy 

resonance, independent of the parameter used to tune it through a threshold. 

This explains previously puzzling results obtained in low-energy reactive scattering. 
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Figure 2. Real (red) and imaginary (green) parts of the scattering length for 3 He + 
NH collisions in the vicinity of an inelastic Feshbach resonance at a kinetic energy of 
1(T 6 K. The lines show the results of Eq.l20l This is the same resonance as shown in 
the right-hand panel of Figure [T] 



Quemener et al. [25] and Cvitas et al. [12J have investigated the sensitivity of scattering 
cross sections in Na + Na2 and Li + Li2 to variations in the potential energy surface. 
Scaling the potential tunes reactive scattering resonances across threshold, and this 
produces oscillations in the elastic and inelastic cross sections. In these systems the 
couplings to individual vibrationally inelastic channels are somewhat reduced by the 
large kinetic energy release, so that for low initial v (with relatively few inelastic 
channels) some significant resonant peaks remain. For initial v = 1, the cross sections 
oscillate by about a factor of 10 as resonances cross threshold. Even this corresponds to 
a relatively small oscillation in the complex scattering length (small a res ). However, 
the amplitudes of the oscillations decrease substantially with increasing vibrational 
excitation of the colliding molecules and are almost smooth for v = 3 for both Na 
+ Na 2 [25j and Li + Li 2 [12]. 

Quite different behaviour has been observed in F + H 2 reactions [26], but is also 
explained by the present theory. Bodo et al. [26] investigated the effect of scaling the 
reduced mass and observed pole-like behaviour in the scattering length and large reactive 
cross sections as a resonance was tuned across threshold. In this case the resonant state 
is localised in the entrance channel of the reaction, while the only exoergic channels 
are reactive ones that are separated from the entrance channel by a high barrier. T^ el 
is thus reduced relative to r^ - Because of this, a res is large (>100 A) and no strong 
suppression of the resonant peaks occurs. 

The considerations of the present paper lead to a remarkable conclusion. It 
has been commonly believed that collision cross sections in the ultracold regime are 
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extremely sensitive to details of the potential energy surface, and that for molecules these 
dependences would be even more limiting than for atoms. The present paper has shown 
that this is true only when the resonant state is coupled much more weakly to inelastic 
(exoergic) channels than to the incoming channel. There are some systems where 
inelastic processes are weak enough for scattering lengths to reach near-infinite values 
at zero-energy resonances. However, in other cases inelastic processes will suppress this 
behaviour. In general terms, the resonant peaks are suppressed by inelastic scattering 
unless there is a specific mechanism that reduces the coupling to inelastic channels. 
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